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We present the first constraints on pure-gravity sector Standard-Model Extension 
(SME) parameters using Lunar Laser Ranging (LLR). LLR measures the round 
trip travel time of light between the Earth and the Moon. With 34-|- years of LLR 
data, we have constrained six independent linear combinations of SME parameters 
at the level of lO''^ to 10"". There is no evidence for Lorentz violation in the 
LLR dataset. 



1. Introduction and Motivation 

Two of us (C.W.S. and J.B.R.B.) are members of the Apache Point Obser- 
vatory Lunar Laser-ranging Operation (APOLLO), a next-generation LLR 
facihty, capable of miUimeter-precision lunar range measurements (see the 
article by T.W. Murphy in these proceedings). The APOLLO project was 
motivated by the realization that an order-of-magnitude improvement in 
fundamental physics constraints (e.g. equivalence principle, gravitomag- 
netism, gravitational law and G, to name a few) could be achieved 
with straightforward improvements to the standard LLR apparatus. 

With the recent description of the pure-gravity sector of the SME,^ we 
learned that LLR can also provide incisive constraints on Lorentz Violation. 
The predicted LLR observable under Lorentz Violation is a periodic per- 
turbation to the Earth-Moon range with the leading order effects occuring 
at four distinct frequencies: 2uj, uj, 2cj — and fi®. Here uj is the lunar or- 
bital (sidereal) frequency, ujq is the anomalistic lunar orbital frequency and 

is the mean Earth orbital (sidereal) frequency. Although the APOLLO 
program is still in the data collection phase, there are more than three 
decades of freely-available archival LLR data on a public archive.^ In this 
article, I present the SME parameter constraints that result from our anal- 
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ysis of archival LLR data. These are the first LLR-based constraints on 
pure-gravity SME parameters. 

2. The LLR Dataset and Analysis Software 

LLR measures the time of flight of photons between a telescope on the Earth 
and corner cubes on the lunar surface."^ LLR data is typically presented in 
"normal points" which are typically generated from a few to a hundred 
lunar signal photons collected over a span of 1 to 5 minutes. Our analysis 
makes use of archival data from September 1969 through December 2003. 

In the analysis of LLR normal points, a set of residual ranges is com- 
puted by subtracting the model's predicted range from the observed range. 
The range sensitivity with respect to each model parameter (the partial 
derivatives) is also computed at the time of each normal point. The resid- 
uals and the partial derivatives are then used to compute optimal model 
parameter values via a weighted linear least-squares fit. For our analy- 
sis, we used the Planetary Ephemeris Program (PEP),^ which is currently 
maintained by one of us (J.F.C.). To our knowledge, it is the only publicly 
available LLR analysis software. 

Typically, the lunar range model is formulated in the parametrized 
post-Newtonian (PPN) framework,^ which permits model-independent con- 
straints on metric theories of gravity. At present, no ephemeris models 
explicitly incorporate SME parameters. You can, however, think of these 
models as implicitly including the SME parameters but with values pegged 
at zero {i.e. no SME perturbation to the lunar orbit). It is therefore only 
necessary to compute, by hand, the partial derivative of range with respect 
to each SME model parameter (see Table 1). The analysis code can then 
provide SME parameter adjustments simultaneously with the other model 
parameters. A covariance matrix including the correlations between the 
SME parameters and all other model parameters is also produced. 

The main drawback to this approach is that one cannot perform an 
iterative analysis in which one takes the best-fit model parameter values 
and uses them to re-integrate the equations of motion to refine the model 
parameter values. We accept this limitation because the non-SME model 
parameters have been highly refined through iterative solutions over the 
past several decades and so the solution sits firmly in the linear regime 
already. Furthermore, the addition of the SME parameters preserves the 
linearity because the lunar range is strictly linear in the SME parameters 
(see Table 2 of Ref. 1), so no iteration is necessary. 
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Table 1. SME parameter partial derivatives. Symbols used here are explained in Ref. 1. 
SME Parameter Partial Derivative of Lunar Range with Respect to SME Parameter 



3. Systematic Errors 

The solar system is complex. Predictions of the lunar range rely on models 
of planetary and asteroid positions, gravitational harmonics of the Sun, 
Earth and Moon and various relativistic and non-gravitational effects (to 
name a few). Solar system models have many hundreds of parameters 
that influence the Earth-Moon range time. There are strong correlations 
between model parameters. As a result, solutions will suffer from systematic 
errors in model parameter estimates that can dominate the formal errors 
reported by the least-squares analysis. In this work, we account for the; 
underestimation of model parameter uncertainties by scaling the formal 
parameter errors reported by the least-squares analysis by a uniform factor, 
F. This is numerically equivalent to uniformly scaling the uncertainty of 
each normal point by F. Essentially, we uniformly down-weight the data. 

The F factor is empirically determined by holding the SME parameter 
values at zero but allowing the PPN values /3 and 7 to vary (see Ref. 4 for 
an explanation of /3 and 7) . We know from existing experiments^ that these 
parameters are consistent with their General Relativity values (/3 = 7 = 1) 
to within a part in 10^ or better. We find that we require F = 20 to ensure 
that we are in accord with these earlier results. 

4. SME Parameter Constraints and Verification 

We constrain the SME parameters under the assumption that General Rela- 
tivity is not violated (we set /? = 7 = 1). The resulting parameter estimates 
and their realistic errors (the formal errors scaled by = 20) are reported 
in Table 2. All SME parameters are within l.hFa of zero. There is no 
evidence for Lorentz violation in the LLR data. The fit quality is shown in 



To verify our implementation of the partial derivatives of lunar range 
with respect to the SME parameters, we generated, by hand, a perturbed 
LLR normal point data set by setting s^^—s^"^ = 9x 10^^*^, a lOFcr deviation 




J cos [(2a) -uJo)t + 26] 
sin [(2w ~t^Q)t + 2ff\ 



Fig. I. 
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Figure 1. The agreement between model and data for the first 34 years of LLR. 

from the best-fit value of this parameter. A fit to this data recovers the 
perturbation: s^^ - 5^2 = [(i + 9) ± 0.9] x 10-^° with the other SME 
parameters unchanged. 

5. Conclusions and Future Prospects 

We have analyzed 34+ years of LLR data and have derived constraints on 
six SME parameters combinations (see Table 2). We find no deviation from 
Lorcntz Symmetry at the 10~^ — 10^^^ levc;!. This work provides the first 
LLR-based constraints of SME parameters. 

There are several ways in which these constraints could be improved. 
First of all, by incorporating auxiliary solar system data (e.g. planetary 
radar ranging) model parameter correlations can be reduced, and F de- 
creased. This would allow for tighter constraints on the SME parameters 
using the same LLR dataset. In addition, APOLLO data, which is about 10 
times more precise than the archival data, will soon be ready for analysis. 
With this improved dataset, we will further tighten the SME parameter 
constraints. 
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Tabic 2. SME parameter estimates 
and their realistic (scaled) uncertain- 
ties {Fa) with F = 20. 



Parameter Estimate 
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